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We report a combined experimental and theoretical study of the spin S — ^ nanomagnet 
Cu5(OH)2(NIPA)4-10H2O (Cus-NIFA). Using thermodynamic, electron spin resonance and ^H nu- 
clear magnetic resonance measurements on one hand, and ab initio density-functional band-structure 
calculations, exact diagonalizations and a strong coupling theory on the other, we derive a micro- 
scopic magnetic model of CU5-NIPA and characterize the spin dynamics of this system. The ele- 
mentary five-fold Cu^''' unit features an hourglass structure of two corner-sharing scalene triangles 
related by inversion symmetry. Our microscopic Heisenberg model comprises one ferromagnetic and 
two antiferromagnetic exchange couplings in each triangle, stabilizing a single spin S = | doublet 
ground state (GS), with an exactly vanishing zero-field splitting (by Kramer's theorem), and a very 
large excitation gap of A ~ 68 K. Thus, Cus-NIFA is a good candidate for achieving long electronic 
spin relaxation (Ti) and coherence {T2) times at low temperatures, in analogy to other nanomag- 
nets with low-spin GS's. Of particular interest is the strongly inhomogeneous distribution of the 
GS magnetic moment over the five Cu^^ spins. This is a purely quantum-mechanical effect since, 
despite the non-frustrated nature of the magnetic couplings, the GS is far from the classical collinear 
ferrimagnetic configuration. Finally, CU5-NIFA is a rare example of a S = | nanomagnet showing 
an enhancement in the nuclear spin-lattice relaxation rate l/Ti at intermediate temperatures. 

PACS numbers: 75.50.Xx, 75.10.Jm, 75.30.Et, 76.60. Jx 



I. INTRODUCTION 

The field of molecular nanomagnets has enjoyed an 
enormous experimental and theoretical activity over 
the last few decades.^ Owing to the nanoscopic size 
of their elementary magnetic units, these compounds 
provide experimental access to a plethora of quantum 
mechanical (QM) effects, including quantum tunneling 
of the magnetization^''^ or the Neel vector,^ quantum 
phase interference,^ level-crossings and magnetization 
plateaux.^ They also allow to probe on the macroscopic 
scale the crossover from quantum to classical physics.^ 
Finally, molecular magnets are promising materials for 
spintronic applications** and quantum computing.^ 

Here we report on the magnetic behavior, mi- 
croscopic magnetic model, and spin dynamics of 
Cu5(OH)2(NIPA)4 • IOH2O, hereinafter referred to as 
CU5-NIPA, where the acronym NIPA stands for the 
5-nitro-isophtalic acid ligand. It is an hourglass-shaped 
molecular magnet comprising five Cu^''" spin-^ ions. The 
ground state (GS) of this magnet has a low spin value 
S=^ with a very large spin gap of A ~ 68 K. Therefore, 
CU5-NIPA behaves as a rigid spin S—^ entity in a wide 
temperature range, and resembles other spin- 2 molecu- 
lar magnets, such as Vq.^'^'^^ Given that both compounds 
comprise s = | spins, we expect similarly long electron 
spin-phonon relaxation times Ti , which allow for the ob- 



servation of rich hysteresis effects in pulsed fields. ^°'^^ 
However, in contrast to Vg, here we do not expect abrupt 
steps in the magnetization curve, ^"^ because the present 
compound features an odd number of half-integer spins, 
thus the zero-field splitting vanishes exactly by Kramer's 
theorem. 

In analogy with other molecular magnets with low- 
spin GS's, such as iron-sulfur clusters, ^'^"^'^ heterometallic 
rings, '^^ iron trimers,^ Vis,^^ and single-molecule mag- 
nets (SMM),^^ we expect that CU5-NIPA manifests also 
long coherence times T2, which is a crucial step towards 
implementations in quantum computing.^'' 

Another attractive feature of CU5-NIPA is the pres- 
ence of two corner-sharing scalene triangles, related 
by inversion symmctry.^^'^^ The spin triangle is the 
most elementary unit for highly frustrated magnetism, ^^ 
while its chirality may induce a finite magnetoelec- 
tric coupling, ^"^'^^ and may also be used as a qubit.^^ 
In molecular magnetism, the spin triangle is com- 
mon in many compounds,^ such as Vis,^^'^^ CU3 
clusters, ^^ the chiral Dya cluster, ^^ the cuboctahedron 
Cui2La8,'^'^ as well as the giant icosidodecahedral ke- 
plerates Mo72Fe3o,^^ W72Fe3o,32 MoyzCrgo,^^ Mo72V3o^'^ 
and W72V3o,^^ which host highly frustrating kagome-like 
physics.'^® 

The distinct feature of the present compound is the 
strong distortion of its regular spin triangles, featuring 
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FIG. 1. (Color online) Left panel: crystal structure of Cus-NIPA showing the stacking of the Cus units (one unit is highlighted 
by gray shading) through the NIPA ligand molecules, into a rigid three-dimensional metal-organic framework, with the shortest 
distance of 8.6 A between the magnetic units. Middle and right panels: structure of the Cus magnetic unit and relevant 
magnetic interactions according to the fit of the experimental susceptibility data (Fig. 13). The numbers 1 — 3 denote the 
crystallographic positions Cul-Cu3, with Cu3 residing at the inversion center. Arrows show the classical ferrimagnetic ground 
state. 



three nonequivalent Cu positions, denoted by Cul, Cu2, 
and Cu3 (see Fig. 1). Each spin triangle is centered by the 
/i3-0H group'^^ providing Cu-O-Cu superexchange path- 
ways, while the carboxyl (C00~) group of the NIPA lig- 
and creates an additional Cul-Cu3 pathway (see Fig. 1, 
middle) . This topology leads to three drastically different 
interactions, one of which is ferromagnetic (FM). Despite 
the fact that these couplings do not compete with each 
other classically, the GS has strong QM character, which 
is partly reflected in a strongly inhomogeneous distribu- 
tion of the magnetic moment over the five Cu^^ spins. 

Of particular interest is our experimental finding of 
an enhancement of the ^H nuclear spin-lattice relaxation 
rate l/Ti at a characteristic temperature slightly below 
the spin gap (T ~ 40 K). While such an enhancement 
has been repeatedly found in numerous antiferromag- 
netic (AFM) homometallic'^^ and heterometallic'^^ rings 
of spins s > ^ , it is very rare for s — ■^. We argue that 
the origin of the peak is the same in both cases, namely 
the slowing down of the phonon-driven magnetization 
fluctuations. '^*'^°"^^ However, there are several qualita- 
tive differences with the case of homometallic rings, re- 
lated to the sparse excitation spectrum of CU5-NIPA and 
the presence of inequivalent Cu^"*" sites. 

The organization of this article is the following. We 
begin in Sec. II with details on the sample prepara- 
tion and methods. Next, we discuss bulk specific heat 
and magnetization (Sec. Ill A), as well as local electron 
spin resonance (Sec. IIIB) and ^H nuclear magnetic res- 
onance data (Sec. IIIC). The microscopic description of 
CU5-NIPA in terms of the isotropic Heisenberg model 
is based on ab initio density-functional band-structure 
calculations (Sec. IV A) and facilitates the evaluation of 
magnetic properties (local magnetizations, nature of the 
GS, etc.) using exact diagonahzations (Sec. IV B) and a 
strong coupling theory (Sec. IV C). Finally, we conclude 
with a discussion and some interesting perspectives of 
this study in Sec. V. 



II. METHODS 

A powder sample of CU5-NIPA was prepared according 
to the procedure described in Ref. 21. Sample purity was 
checked by powder x-ray diffraction (Huber G670 Guinier 
Camera, CuKcji radiation, 26* = 3— 100° angular range). 

The magnetic susceptibility (x) was measured in the 
temperature range 1.8 K < T < 400 K in an applied field 
of fioH = 1 T using the commercial Quantum Design 
MPMS SQUID. The magnetization isotherm M vs. H 
was measured at T = 2 K in fields up to 14 T using the 
vibrating sample magnetometer (VSM) option of Quan- 
tum Design PPMS. Additionally, pulsed-field measure- 
ments in fields up to 60 T were performed at 1.4 K in the 
Dresden High Magnetic Field Laboratory. Details of the 
measurement procedure are described in Ref. 43. 

The heat capacity Cp(T) was measured on a 
small piece of pellet over the temperature range 
2 K < r < 300 K in zero field using the Quantum De- 
sign PPMS. 

The electron spin resonance (ESR) measurements were 
performed at Q-band frequencies (/ — 34 GHz) using a 
standard spectrometer together with a He-flow cryostat 
that allows us to vary the temperature from 1.6 to 300 K. 
ESR probes the absorbed power P of a transversal mag- 
netic microwave field as a function of a static and exter- 
nal magnetic field B. To improve the signal-to- noise ra- 
tio, we used a lock-in technique by modulating the static 
field, which yields the derivative of the resonance signal 
dP/dB. 

The nuclear magnetic resonance (NMR) measurements 
were carried out using pulsed NMR technique on ^H nu- 
clei (nuclear spin / = 5 and gyromagnetic ratio 7„/27r = 
42.576 MHz/T) in the 2-230 K temperature range. The 
NMR measurements were performed at two different ra- 
dio frequencies of 70 MHz and 38.5 MHz, which corre- 
spond to an applied field of about 1.6608 T and 0.9135 T, 
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FIG. 2. (Color online) Left panel: Inverse magnetic susceptibility of Cus-NIPA measured upon heating from 2 K (circles, 
hydrated form) and upon cooling from 400 K (triangles, dehydrated form). Lines show Curie- Weiss fits to Eq. (1), as described 



in the text. Right panel: x'T plot showing the formation of the j|) state in the pristine sample {x'T 



and the high- 



temperature paramagnetic regime of the dehydrated sample (x*T = |). Here, X* ~ X 
according to ESR (Sec. IIIB). For the x{T) plot, see Fig. 13. 
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N = Na/^, and g = 2.2 



respectively. Spectra were obtained by Fourier transform 
of the NMR echo signal. The ^H spin-lattice relaxation 
rate, 1/Ti, was measured by using the conventional sat- 
uration pulse sequence. 

The electronic structure of CU5-NIPA was calculated in 
the framework of density functional theory (DFT) using 
the VASP code'^'* for crystal structure optimization and 
FPLD'*^ for the evaluation of magnetic parameters. The 
generalized gradient approximation (GGA)^^ exchange- 
correlation potential was augmented with the mean-field 
DFT-|-[/ correction for Coulomb correlations in the Cu 
2>d shell. The DFT+f7 parameters were chosen as Ud = 
9.5 eV (on-site Coulomb repulsion), Jd = 1 eV (on-site 
Hund's exchange), and fully- localized-limit (FLL) fla- 
vor of the double-counting correction, following earlier 
studies of Cu^^-based magnets.'*^ Reference calculations 
with other exchange-correlation potentials and double- 
counting corrections arrived at qualitatively similar re- 
sults. The reciprocal space was sampled by a fc mesh 
with 64 points in the first Brillouin zone. The conver- 
gence with respect to the k mesh was carefully checked. 
Details of the computational procedure are reported in 
Sec. IV A. 

Thermodynamic and GS properties of the Cus 
molecule were evaluated by full numerical diagonaliza- 
tions. 



III. EXPERIMENTAL RESULTS 
A. Thermodynamic properties 

1. Magnetization 

The magnetic susceptibility of CU5-NIPA was mea- 
sured for the as-prepared sample under field-cooling (FC) 
condition upon heating from 1.8 K to 400 K. At 400 K, 



the sample was kept inside the SQUID magnetometer for 
about 20 minutes, and xC^) '^^^ measured again upon 
cooling from 400 K to 1.8 K. The drastic difference be- 
tween the heating and cooling curves (Fig. 2) indicates 
the decomposition of CU5-NIPA shortly above room tem- 
perature. Indeed, the sample color changed from blue to 
green, and the weight loss of 13% was detected. This 
weight loss is in good agreement with the preceding ther- 
mogravimetric data of Ref. 22 that reported the weight 
loss of 13.6% at 380 - 400 K. Although the decomposi- 
tion process is tentatively ascribed to the release of water 
molecules (expected weight loss of 13.15 %)^^, no detailed 
information on the decomposed sample is available in the 
literature. 

The pristine sample shows a sharp decrease in the mag- 
netic susceptibility (increase in 1/x) upon heating and a 
bend around 30 K followed by the nearly linear regime of 
1/x above 200 — 250 K. In contrast, the decomposed sam- 
ple remains paramagnetic down to at least 10 K. Below 
320 K, where the decomposition process is finished, the 
inverse susceptibility of dehydrated CU5-NIPA follows a 
straight line that can be fitted with the Curie- Weiss law 
in the 20 — 320 K temperature range, ^* 
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The resulting Curie constant C ~ 0.383 emu K/(mol Cu) 
leads to an effective moment fi^s — 1-75 /is, which is 
close to 1.73 /is expected for spin-i. The Weiss temper- 
ature is 6* ~ 9 K. 

The inverse susceptibility of the pristine CU5-NIPA 
sample does not have a well-defined paramagnetic re- 
gion, because the data above 320 — 350 K are affected 
by the decomposition. A tentative Curie- Weiss fit in 
the 220 — 320 K range yields an effective moment of 
fics — 1-75 hb (C = 0.383 emu K/(mol Cu)), which is 
same as in the dehydrated sample. However, the Weiss 
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FIG. 3. (Color online) Magnetization curve of Cus-NIPA 
measured at 1.8 K up to 14 T. The dashed line denotes the 
saturation of the ||) state at Ms = 1.19 ^b/^-Vl. The solid 
line is the Brillouin function calculated at 1.8 K with g — 2.38. 
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FIG. 4. (Color online) Specific heat of Cus-NIPA divided 
by temperature (Cp/T). The inset in the upper left corner 
magnifies the data at low temperatures. The inset in the 
bottom right corner shows the specific heat (Cp). 



temperature of ~ 34 K in the pristine sample is notably 
larger than 6* ~ 9 K observed in the dehydrated sample. 

The conspicuous difference between the hydrated and 
dehydrated samples suggests a huge effect of dehydra- 
tion on the magnetism of the Cu5 molecule. The dehy- 
dration proceeds in a single step and results in the loss 
of all 10 water molecules per formula unit. As 4 out 
of these 10 molecules enter the first coordination sphere 
of Cu (Fig. 1, left), a large effect on the local environ- 
ment of Cu sites and, thus, on the magnetism should 
be expected. ^^ Unfortunately, no structural data for the 
dehydrated version of Cus-NIPA are presently available. 
Therefore, we restrict ourselves to a detailed study of the 
hydrated compound. 

At low temperatures, the inverse susceptibility of Cus- 
NIPA approaches another linear regime, with the Curie 
constant C = 0.101 emu K/(mol Cu) (/Xeff ~ 0.90 ^s) 
and vanishingly small 6 ~ 0.5 K. This low-temperature 
paramagnetic state can be also seen in the x*T plot 
(x* — X (-^ff^Ms/^s) is the reduced susceptibihty), 
where x*T approaches the value of \ expected for the 
total spin -S* = ^ per molecule (Fig. 2, right). Likewise, 
the magnetization curve of Cus-NIPA (Fig. 3) saturates 
at Ms = 1.19 yitB /molecule, which is, however, higher 
than the free-electron value of 1 /x^/f.u. This difference 
can be well accounted for by the large ^-value of 2.38 ac- 
cording to Ms = gSfiB- The same g- value explains the 
low-temperature effective moment: C = Ng^fi^g/'dkB = 
0.106 emu K/(mol Cu), where we use N = iVyi/5 for 
the magnetic moment of ^ per Cus molecule. As we ex- 
plain later, this change in the 5- value (Sec. IV B) can be 
traced back to the presence of three non-equivalent Cu 
positions with different (7-tensor anisotropics. At room 
temperature, the Cu spins are nearly independent, and 
their g-values average with same weights for all Cu posi- 
tions. At low temperatures, the powder averaging of the 
(7- values is determined by the distribution of the magne- 



Our results suggest that at low temperatures Cus- 
NIPA is in the paramagnetic state with the total moment 
of S* = i per Cu5 molecule. This state is further denoted 
as \\) ■ The experimental magnetization curve follows 
the Brillouin function 



B{H) = g^ieS x th 






(2) 



tization in the S 



GS. 



with S' = i, 5 = 2.38, and T = 1.8 K (Fig. 3). A marginal 
departure of the experimental curve toward higher fields 
may be due to anisotropics and/or intermolecular cou- 
plings. The maximum deviation of about 0.3 T puts an 
upper limit of about 0.5 K on the total energy of such 
couplings. This low energy scale is in agreement with the 
large spatial separation between the molecules (Fig. 1, 
right). 

Below 4 K, the susceptibility of CU5-NIPA departs 
from x*T = \ (Fig. 2, right). This decrease in x*T may 
indicate an evolution of the system toward a long-range 
magnetic order between the CU5 molecules. However, we 
were unable to see any clear signatures of a magnetic 
transition in the susceptibility data measured down to 
2 K. 



2. Specific heat 

The specific heat (Cp) of CU5-NIPA is smooth down 
to 2 K (Fig. 4). It steeply increases up to 60 — 70 K 
and shows a slower increase at higher temperatures, as 
shown by a broad maximum in the temperature depen- 
dence of Cp/T . The signal is dominated by the phonon 
contribution, whereas the magnetic part is quite small. 
The total magnetic entropy of CU5-NIPA is 5i?ln2 ~ 
28.8 J mol"'^ K~^, which is only 2 % of the total entropy 
of about 1100 J mor^ K'^ released up to 200 K (the 
latter is obtained by integrating the temperature depen- 
dence of Cp/T). 



At 200 K, the heat capacity of CU5-NIPA is still very 
far from its maximal Dulong-Petit value of Cp = 3RN ~ 
2767 J mol^^ K~^, because CU5-NIPA features a large 
number of high-energy phonon modes related to the 0-H, 
C-C, C-N, and C-0 vibrations. The entropy associated 
with these modes can be released at very high temper- 
atures, only. The complexity of the CU5-NIPA struc- 
ture and the presence of multiple phonon modes of differ- 
ent nature hinder a quantitative analysis of the specific 
heat data. A non-magnetic reference compound would 
be ideal to extract the magnetic contribution and ana- 
lyze it in more detail. Unfortunately, such a reference 
compound is presently not available. 

At low temperatures, the heat capacity of CU5-NIPA 
does not fall smoothly to zero. The temperature depen- 
dence of Cp/T shows an increase below 3 K (see the up- 
per left inset of Fig. 2). The origin of this behavior is 
presently unclear. The increase in Cp/T could signify a 
Schottky anomaly or a proximity to the magnetic order- 
ing transition. However, our susceptibility data, as well 
as NMR (Sec. IIIC), rule out any magnetic transition 
down to 2 K. The heat capacity data are consistent with 
these observations. 



B. ESR 

The ESR spectrum of CU5-NIPA could be detected 
from the lowest measured temperature (5 K) up to 135 K. 
Above 135 K, the ESR line broadens and eventually be- 
comes invisible. Between 70 K and 135 K, the spec- 
tra contain only one line that can be fitted with a sin- 
gle Lorentzian function (IL-Fit) providing the ESR pa- 
rameters, the linewidth AS and q- factor, q — — '^ — , 
where iJros is the resonance field. Below 70 K, the sin- 
gle ESR line splits, so that at 40 — 60 K the fit is only 
possible with two Lorentzian lines (2L-Fit), whereas at 
even lower temperatures a plethora of narrow lines ap- 
pear (Fig. 5). At 5 K, more than 20 lines are visible 
in the spectrum. This complex structure is typical in 
low-temperature ESR spectra of molecular magnets and 
is related to the hyperfine coupling. ^"'^^ The complexity 
of the signal and the presence of multiple Cu sites and 
anisotropy parameters prevent us from a detailed anal- 
ysis of the low-temperature powder spectra. One would 
need much higher frequencies and, preferably, single crys- 
tals of CU5-NIPA, in order to discriminate the multiple 
resonances that are visible below 40 K. So we restrict 
ourselves to the analysis of the data above 40 K. 

Fig. 6 shows the temperature evolution of the ESR g- 
factor, linewidth, and line intensities. Above 70 K, the 
observed powder-averaged values of 5 ~ 2.2 are typical 
for Cu^+ in the planar oxygen environment. ^^'^"'^^ Our 
microscopic insight into the energy spectrum of the CU5 
molecule (Sec. IV B) suggests that the formation of two 
lines below 70 K is related to two lowest S — ^ energy 
levels, which are separated by A ~ 68 K. Indeed, lines 
1 and 2 show different temperature evolution. While the 
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FIG. 5. (Color online) Typical ESR signal of Cus-NIPA at 
5 K (solid line) and 65 K (dashed line). For better visibility, 
the intensity of the high-temperature signal was increased by 
a factor of 20. 
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FIG. 6. (Color online) Temperature dependence of the ESR 
signals: g- factor (top), linewidth (middle), intensity (bot- 
tom). The data at higher temperatures (diamonds) were ob- 
tained by fitting the spectra with one Lorentzian line (IL-Fit) . 
Open circles and triangles show the data from fitting the spec- 
tra with two Lorentzian lines (2L-Fit). Representative error 
bars are shown for several data points, only. All intensities 
are scaled to the intensity of line 2 at 65 K. 



intensity of line 2 is reduced upon cooling, the intensity 
of line 1 is growing. In Sec. IV B below, we provide the 
explicit dependence of the effective g-tensors for the GS 
and the first excited doublet (in terms of the individual 
g-tensors of the three inequivalent Cu sites) in order to 
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FIG. 7. (Color online) Top panel: Fourier-transform ^H 
NMR spectra measured at different temperatures. Arrows 
show the shifts of the left and right shoulders upon cooling. 
Note that while the right shoulder is moving toward higher 
frequencies, the left shoulder shows a non-monotonic behav- 
ior. Bottom panel: Full width at half-maximum (FWHM) 
of ^H NMR spectra plotted as a function of temperature (T) 
measured in an external field hqH — 1.6608 T. 



demonstrate the origin of their difference. 



C. ^H NMR 



Linewidth 



Considering the nuclear spin / = ^ of the ^H nu- 
cleus, one expects a single spectral line for each of the 
17 nonequivalent proton sites (see Table III of App. B). 
However, these lines merge into a single broad line with 
a nearly Gaussian line shape over the whole measured 
temperature range (Fig. 7, top). The linewidth at 215 K 
is 95 kHz and almost comparable with the linewidth re- 
ported for other molecular magnets (see, e.g., Ref. 52). 
The line position shifts weakly with temperature. The 
bottom panel of Fig. 7 shows the temperature evolu- 
tion of the full width at half-maximum (FWHM) of the 
^H NMR line measured in an applied field oi ^qH = 
1.6608 T. It is almost temperature-independent at high 
temperatures and increases progressively upon cooling 
below about 30 K. 

The shape and width of the ^H NMR spectra are gov- 
erned by two main interactions: the nuclear- nuclear dipo- 
lar interactions and the hyperfine couplings between the 
proton and unpaired electrons at the Cu sites. Therefore, 



we can write FWHM as:^^ ^^ 

FWHM ex ^/{Aiy^)d + {Aiy^) 



(3) 



where the broadening due to the nuclear-nuclear dipolar 
interactions ((Az/^)^) is temperature-independent, while 
the broadening due to the hyperfine couplings ((Ai/^)„i) 
scales with the local susceptibilities. For each given pro- 
ton site p. 
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where A- is the dipolar coupling constant between the 
proton and the Cu^+ ions at site j, and Xj is the lo- 
cal susceptibility. We should note here that, in general, 
the temperature dependence of Xj is different from that 
of the bulk susceptibility x (Fig- 2), especially at higher 
temperatures, and they are also different from each other 
(see Fig. 15). The sharp increase in FWHM at low 
temperatures can be ascribed to the corresponding low- 
temperature increase in Xj- 

The shallow minimum in FWHM observed around 
70 K is paralleled by the peculiar evolution of the spec- 
tral line (Fig. 7, top) . On cooling down, the right side of 
the line shifts weakly but monotonously to the right, fol- 
lowing the central position of the line. The left side, on 
the other hand, shows a non-monotonic behavior, shifting 
backwards for the data at 20 K and below. The origin of 
this feature (and the minimum in FWHM) can in princi- 
ple be traced back to the temperature dependence of the 
local Cu^+ moments, and indeed such a non-monotonic 
behavior is shown by the moments on the Cu2 sites. As 
we discuss in detail below (Sec. IV B and Fig. 15), at 
T ^ A ~ 68 K, the moments of the two Cu2 sites are 
antiparallel to the field, owing to the large negative ex- 
change field exerted by the neighboring Cul and Cu3 
sites. These exchange fields are balanced by entropy at 
some characteristic temperature T*, at which the Cu2 
moments turn positive. Our calculations based on the 
actual exchange couplings give T* ~ 38 K (see Fig. 15 
below). At even higher temperatures, the Cu2 moments 
attain their paramagnetic Curie-like behavior of isolated 
spins. The contribution to the second moment from the 
Cu2 sites is then expected to decrease down to zero and 
then increase again as we cool down, starting from the 
high-temperature paramagnetic to the low-temperature 
"ferrimagnetic" S = ^ state of Cus-NIPA. 



2. Wipe-out effect 

The wipe-out effect is common in molecular 
nanomagnets,'^^'^^"^^ and refers to the gradual loss 
in the NMR signal intensity below a characteristic 
temperature. Here, we have measured the temperature 
dependence of Mxy{0)T, where Mxy{0) is the transverse 
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FIG. 8. (Color online) Normalized integrated intensity 

(/r//niaxrniax) measured at hqH — 0.9135 T. The strong 
drop in the intensity below 100 K is the wipe-out effect. The 
line shows the maximum intensity (/maxTmax) attained at high 
temperatures. 



magnetization at time t = 0, obtained by the extrapola- 
tion of the Mxy{t) recovery curve to t = 0. The integrated 
intensity Mxy{0)T is proportional to the total number 
of protons resonating at the irradiated frequency. The 
normalized integrated intensity of the NMR signal in 
CU5-NIPA as a function of temperature is shown in 
Fig. 8. We find that the loss of the NMR signal begins 
below 150 K and becomes more pronounced at low 
temperatures. The onset of the NMR signal loss around 
150 K coincides with the regime where l/Ti starts 
increasing towards the maximum (see Fig. 9 below). 

The loss in the NMR signal is related to the slowing 
down of the electron spin dynamics which, as we explain 
below, is also responsible for the enhancement in 1/Ti. 
On decreasing T, a progressively larger fraction of pro- 
tons close to the magnetic Cu^^ ions attain spin-spin 
relaxation times (T2) shorter than the dead time (rd) of 
the spectrometer, hence the signal from these protons 
cannot be detected. At low temperatures, only the pro- 
tons that are very far from the magnetic ions contribute 
to the signal intensity. For further details, see Refs. 56 
and 57. 



3. Nuclear spin-lattice relaxation rate 1/Ti 

The spin- lattice relaxation rate 1/Ti for CU5-NIPA was 
measured at two different applied fields. The recovery of 
the longitudinal nuclear magnetization after a saturation 
pulse was fitted well by the stretched exponential func- 
tion (top panel of Fig. 9, inset) 
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where M{t) is the nuclear magnetization at a time t after 
the saturation pulse, and Mq is the equilibrium magne- 
tization. The value of the exponent /3 for both fields was 
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FIG. 9. (Color online) Top panel: Spin-lattice relaxation 
rate 1/Ti vs. T measured at two different values of external 
field, 0.9135 T and 1.6608 T. The inset shows the longitudi- 
nal recovery curves at four representative temperatures; solid 
lines are the fits using Eq. (5). Bottom panel: l/{xTiT) 
vs. T for two different magnetic fields. The inset shows the 
temperature-dependent ujc (decay rate of the total moment 
Sz) extracted from the 1/Ti data; the solid line is the fit 
using Eq. (7). 



found to decrease slowly from 0.95 to 0.7 upon cooling. 
The deviation of /3 from unity refiects the distribution of 
relaxation rates according to different hyperfine couplings 
between the ^H nuclei and Cu^"*" ions.^^'^^'^" Accordingly, 
the spin-lattice relaxation rate l/Ti obtained by fitting 
the recovery with Eq. (5) provides a value averaged over 
the whole distribution. 

The T dependence of l/Ti is presented in the top panel 
of Fig. 9 for two values of the external field. At high tem- 
peratures (T > 150 K), 1/Ti is almost T-independent. 
This behavior is typical for uncorrelated paramagnetic 
moments fiuctuating fast and at random. ^^'^^ At lower 
temperatures, 1/Ti increases and passes through a max- 
imum at r ~ 40 K. Such a characteristic enhancement 
has been found in numerous AFM homometallic'^* and 
heterometallic'^^ rings built of spins s > 5 , but in spin- ^ 
systems it is very rare. The only example known to us 
is the Cug magnet with a high-spin 5 = 3 GS.®'^'®^ Sys- 
tems with predominantly AFM couplings and low-spin 
GS (e.g., V12 having an 5 = GS, Ref. 65) do not show 
such a feature in l/Ti. 

In homometallic rings, the maximum in 1/Ti essen- 
tially signals the slowing down of phonon-driven spin 
fiuctuations.'^^'^^"'*^'^^ The equivalence between the spins 
(by virtue of the nearly perfect rotation symmetry of the 
ring) allows to express 1/Ti in terms of the spectral den- 
sity of the total magnetic moment 5*^ of the molecule."'" 
As shown numerically^'^ and by a microscopic theory,^^'^^ 



the phonon-driven decay of 5*^ proceeds independently 
from the remaining observables of the problem, which in 
turn leads to a single Lorentzian form for l/Ti: 
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where A is the average square of the transverse hyperfine 
field, ujl is the nuclear Larmor frequency, and uj^ = -;^rhps 
is the decay rate of the total moment Sz ■ The enhance- 
ment in the spin-lattice relaxation rate 1/Ti takes place 
when uJc approaches the order of magnitude of uj^ ■ 

There are two qualitative differences between Cus- 
NIPA and homometallic rings that should be emphasized, 
though. First, Cus-NIPA features a very sparse excita- 
tion spectrum (Sec. IV B), which means that there are 
very few spin-phonon channels available for relaxation. 
According to Ref. 42, this also means that here Wc is 
not expected to show the strong power-law ^ T^ or T* 
behavior. Instead, one expects a much weaker T depen- 
dence, with longer relaxational times in a wide tempera- 
ture range. 

The second difference is the fact that the Cus molecule 
comprises three inequivalent Cu sites with different local 
magnetizations, and so the above single-Lorentzian for- 
mula for 1/Ti is not valid any longer, but instead a multi- 
Lorentzian form should be expected on general grounds, 
as in the case of the heterometallic CryNi ring.^^ 

Given the large spin gap, A ~ 68 K, one could still 
argue in favor of the single-Lorentzian formula of Eq. (G) 
over a wide low-T range, since the Wigner-Eckart theo- 
rem allows to replace individual spin operators with the 
total spin operators times a constant. However, such a 
treatment would only capture the resonant spin-phonon 
transitions between the two Zeeman-split levels of the 
GS, but these transitions are prohibited by Kramer's the- 
orem. Instead, for temperatures well below the spin gap 
A ~ 68 K, the relaxation dynamics of the spins must 
be controlled by inter-multiplet Orbach transitions^^ be- 
tween the GS doublet and the lowest magnetic excitation 
with the energy A (see below). Again, this is quite sim- 
ilar to CryNi, with the only difference that the spin gap 
is much smaller there (~ 14 K).^^ 

With the above remarks in mind, we may still use the 
above single-Lorentzian formula for 1/Ti to extract the 
temperature dependence of uj^, but the latter is now a 
representative measure of the relaxation dynamics (as 
probed by NMR), and not necessarily the relaxation rate 
of the total moment. For simplicity, we have plotted 
1/xTiT vs. T in the bottom panel of Fig. 9 that clearly 
shows a broad maximum reflecting the slowing down of 
fluctuating moments. Following Eq. (6) and the pro- 
cedure outlined in Ref. 38, we extracted lOc from the 
l/xTiT data at ^qH = 0.9135 T. The resulting T de- 
pendence of UJc (shown in the lower inset of Fig. 9) con- 
firms that it is much weaker compared to the strong 
power-law dependence found in homometallic rings. "^^ In 
particular, there is a wide low-T range over which uJc 
follows the inter-multiplet Orbach relaxation processes 




Energy (eV) 

FIG. 10. (Color online) LDA density of states for Cug-NIPA. 
The shading shows the contribution of Cu orbitals. The inset 
magnifies five bands at the Fermi level {E = 0), with the 
atomic contributions denoted by the dashed (Cu) and solid 
(O) fines. 



mentioned above. Following Ref. 67, 
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where pm = 2043 kg/m^ is the mass density,^^ v is the 
sound velocity in CU5-NIPA (typical value c ~ 1500 
ni/scc in nanomagnets^^), and A stands for the spin- 
phonon coupling energy parameter, which is related e.g. 
to the fluctuating portion of the Dzyalozinskii-Moriya in- 
teractions present in this system. Note that we neglect 
the Zeeman splitting contribution to the resonance en- 
ergy, because it is negligible compared to A. By fltting 
the Wc data with Eq. (7), we obtain X/ks — 3.9(ti/c)2 K. 
A more complete quantitative understanding of 1/Ti 
(e.g. the high-temperature behavior and held depen- 
dence) must take into account the multi-exponential be- 
havior of the relaxation, discussed above, but also the 
presence of the wipe-out effect discussed above, see e.g. 
Ref. 66. 



IV. THEORY 

A. Microscopic magnetic model 

To determine individual exchange parameters in the 
Cu5 molecule, we calculate the electronic structure of 
the CU5-NIPA compound. This procedure requires reli- 
able crystallographic information, including precise posi- 
tions of all atoms in the unit cell. However, the available 
structural data^^'^^ are obtained from x-ray diffraction 
that has only limited sensitivity to the positions of hy- 
drogen atoms. Therefore, we used the literature data as 
a starting model and optimized the hydrogen positions, 
whereas all other atoms were kept flxed. The equilibrium 
positions of hydrogen are listed in Table HI of App. B. 
The relaxed structure is 18.6 eV/f.u. lower in energy 




FIG. 11. (Color online) LDA band structure of Cus- 

NIPA (thin light lines) and the fit with the tight-binding 
model (thick dark lines). The Fermi level is at zero en- 
ergy (white dashed line). The notation of k points is as fol- 

lows: r(o,o,o), x(i,o,o), M(i,i,o), y(o,i,o), ^(0,0,1), 

r(|,0, I), -R(|, |, I), and A{0, |, |), where the coordinates 
are given in units of the reciprocal lattice parameters. 



TABLE I. Interatomic distances d (in A), Cu-O-Cu bridging 
angles ip (in deg), hopping parameters ti (in meV), and ex- 
change integrals Ji (in K) in Cus-NIPA. The Ji values are 
obtained from DFT+U calculations. The AFM contribu- 
tions Ji ^^ are evaluated as Atl/UeS, where [/eft is the ef- 
fective on-site Coulomb repulsion. The FM contributions are 

FM _ J _ tAFM 
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than the starting model taken from the literature. This 
energy reduction should be ascribed to the elongation of 
0-H distances that are unrealistically short (about 0.8 A) 
in the experimental structural data.^^ 

The LDA energy spectrum of CU5-NIPA (Fig. 10) com- 
prises narrow lines that represent molecular orbitals of 
the NIPA molecules and the Cu5(OH)2 unit. The states 
at the Fermi level belong to the narrow Cu 3d bands with 
a sizable admixture of O 2p. The local environment of Cu 
atoms resembles the conventional CUO4 plaquette units 
(four- fold coordination, see Fig. 1, left). Therefore, the 
highest-lying Cu 3d bands have predominantly x'^ — y^ 
origin, following the crystal-field levels of Cu^"'" with x 
and y axes lying in the plane of the CUO4 plaquette. 
The five Cu atoms of the CU5 molecule (one molecule 
per unit cell) give rise to five bands, with the middle 
band crossing the Fermi level (Figs. 10 and 11). This 
spurious metallicity is due to the strong underestimate 
of electronic correlations in LDA. DFT-|-f7 calculations 
reveal the robust insulating behavior with a band gap of 
about 2.2 eV in reasonable agreement with the light-blue 
color of the sample. 




FIG. 12. (Color online) Cu d^2_y2-based Wannier function 
for the Cul site in CU5-NIPA. Note that the ligands (COO", 
/i3-OH~, H2O) affect the positions of oxygen atoms and their 
orbitals. However, only the orbitals of four oxygen atoms 
contribute to the Wannier function. 



To evaluate the magnetic couplings, we fit the five Cu 
bands with a tight-binding model (Fig. 11) and extract 
the relevant hopping parameters ti using Wannier func- 
tions (WFs) based on the Cu d^2_y2 orbital character. ^^ 
The ti's arc further introduced into an effective Hub- 
bard model with the on-site Coulomb repulsion Ucs- As 
the conditions of the half-filling and strong correlations 
{ti <C Ucff) are fulfilled, the lowest-lying excitations can 
be described by a Heisenberg model with the AFM ex- 
change J,^™ = Atj/U^Q. The values of ti and J^^^ 
obtained from the LDA band structure are listed in Ta- 
ble I and provide an overview of the magnetic couplings 
in CU5-NIPA. Considering only nearest-neighbor interac- 
tions, we find a strong coupling between Cu2 and Cu3, a 
somewhat weaker coupling between Cul and Cu2, and a 
relatively weak coupling between Cul and Cu3. Long- 
range couplings are several times smaller than J^^^. 
The long-range exchanges in the Cus molecule are within 
10 K, whereas the interactions between the molecules are 
below 3 K (note the experimental estimate of 0.5 K in 
Sec. Ill A). As the long-range couplings are much weaker 
than J12, Ji3, and J23, we further restrict ourselves to 
the minimum microscopic model that comprises three 
nearest-neighbor interactions, only. 

Possible FM contributions to the short-range cou- 
plings require that the evaluation of jf"^^ is supplied 
by an independent estimate of total exchange couplings 
Ji. In a so-called supercell approach, total energies of 
coUinear spin configurations are mapped onto the Heisen- 
berg model to yield the Ji values, as listed in Table I. 
All nearest-neighbor interactions have sizable FM compo- 
nents Jf^ — Ji~jf^^^ that reduce J12 and J23, whereas 
J13 eventually becomes ferromagnetic. This way, the tri- 
angular units in the CU5 molecule feature a combination 
of FM interaction J13 and AFM interactions J12 and J23 
(Fig. 1, left). 

The microscopic origin of the magnetic couplings can 
be understood from the analysis of Cu-based Wannier 
functions (Fig. 12). Each WF features the Cu dx'2-y2 
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FIG. 13. (Color online) Fit of the magnetic suscepti- 

bility (x) and inverse magnetic susceptibility (x~^, inset) 
with the exact-diagonalization result for the Cus pentamer 
(Ji2 = — Ji3 = 65 K and J23/J12 = 3.5). 



orbitals together with the 2p orbitals of the surround- 
ing oxygen atoms. These oxygen atoms have different 
chemical environment and belong to one of the ligands: 
^3-OH-, COO", and H2O. While the ligands affect the 
positions of oxygen atoms (note the slight downward dis- 
placement of the O atoms in C00~ groups) and modify 
the "shape" of the p orbitals (note the oxygen atom be- 
longing to the H2O molecule), they do not provide any 
sizable long-range contributions to the WFs. The short- 
range nature of the WFs underlies the diminutively small 
long-range exchange in CU5-NIPA. Further microscopic 
aspects of the magnetic interactions in CU5-NIPA are dis- 
cussed in Sec. V. 

Our microscopic model can be directly compared to the 
experimental data. Using full diagonalization for the Cus 
molecule, we calculate thermodynamic properties and re- 
fine the Ji parameters as to match the experiment. This 
way, we fit the experimental magnetic susceptibility curve 
with J12 = -Ji3 = 62 K and J23 = 217 K (so that 
J12/J23 — f ), as well a.s g = 2.22 and the temperature- 
independent contribution xo = — 1.7 x 10""' emu/mol Cu 
(Fig. 13). The experimental estimates of Ji are in excel- 
lent agreement with the DFT results (Table I), whereas 
the g-value matches the high-temperature g ~ 2.2 from 
ESR (Sec. IIIB). Regarding the sizable Xoi it may origi- 
nate from the diamagnetic gelatin capsule that was used 
in the susceptibility measurement. 

The combination of the FM coupling J13 and AFM 
couplings J12 and J23 implies that the spin triangles in 
Cu5 are non-frustrated. In a classical picture, all three 
couplings are satisfied in a configuration with Cu3 and 
Cul spins pointing up and Cu2 spins pointing down, thus 
leading to the total moment of | per molecule. Indeed, 
at low temperatures CU5-NIPA shows paramagnetic be- 
havior with the magnetic moment of ^ per molecule ( ^ ) 
state). In high magnetic fields, the ||) and, eventually, 
1 1 ) states will be stabilized. The formation of states with 
the higher magnetic moment can be seen from the sim- 
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FIG. 14. (Color online) Upper panel: magnetization of 
CU5-NIPA measured in pulsed magnetic field (in arb. units). 
Bottom panel: magnetization measured in static field up to 
14 T and the simulated curve for J12 = — J13 = 65 K and 
J23/J12 = 3.5 (both in units of ^s/f.u.) The arrow denotes 
the bend in the pulsed-field curve that matches the transition 
between the 1 1 ) and 1 1 ) states of the Cus molecule. 



ulated magnetization curve of CU5-NIPA. The I ^ ) state 
persists up to about 50 T, where the magnetization in- 
creases abruptly until the 1 1 ) state is reached. 

We endeavored to verify this prediction with magneti- 
zation measurements in pulsed fields. Unfortunately, the 
results are strongly affected by dynamic effects. The typ- 
ical duration of the pulse (about 20 /is) is insufficient to 
change the magnetic moment of the molecule. Therefore, 
all features are blurred and, moreover, the curve shows 
a non-zero slope above 5 T (Fig. 14, top), in contrast to 
the flat plateau that is observed in static field (Fig. 14, 
bottom). The apparent mismatch between the pulsed- 
field and static-field data also prevents us from scaling 
the high-field curve. Nevertheless, the bend observed at 
50 — 55 T is likely the signature of the | ^ ) — > 1 1 ) tran- 
sition and confirms our expectations. 



B. Exact Diagonalization 

We are now ready to evaluate the magnetic energy 
spectrum of each Cus molecule and deduce its basic prop- 
erties, according to the DFT values of the exchange cou- 
plings Ji. Using a site labeling convention in accord with 
the site symmetries, the spin Hamiltonian (disregarding 
anisotropics) reads: 

"H^ J23S3-S22' + J12 (81-82 + Si/ -82' — 83-811') , (8) 

where 822' = 82 + 82' , 811' = 8i-f8i', and we have used 
the relation J13 — — J12. In addition to the full SU(2) 
spin rotation group, this Hamiltonian is also invariant 
under the inversion operation (in real space) through the 
central site 83, which maps 83 O 83, 81 o 81', and 
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S2 <->■ S2' . So all eigenstates can be characterized by the 
total spin S, its projection M along some quantization 
axis z, and the parity p (even or odd) under inversion. 

Global spectral structure. A straightforward numerical 
diagonalization of "H yields the energy spectrum given in 
Table II. For each level we also provide the good quantum 
numbers S and p, as well as the expectation values of the 
square of the composite spins S22' , S322' = S3+S22' , and 
Sii'. The essential finding is that the GS is a total spin 
S — ^ doublet and the lowest excitation (also a S = ^ 
doublet) has a very high energy gap A = 0.313J23~68 K. 
Thus the Cus molecule behaves as a rigid S = ^ entity 
for T ^ A, confirming the previous experimental picture 
from magnetization measurements (Sec. IIIAl). 

Looking at Table II, we also find that the spectrum is 
organized in three compact groups of states, which are 
separated by horizontal lines. All states belonging to a 
given group show a very similar value of (S322/), suggest- 
ing that the trimer of S3, S2, and S2' plays a special 
role in the physics of CU5-NIPA. Moreover, the expecta- 
tion values of 832/, S322/, and Sf^, are for all states very 
close to the formula S{S + 1) for some integer or half- 
integer S, suggesting that these composite spins are al- 
most conserved quantities. Below, we shall demonstrate 
that all these spectral features as well as the GS proper- 
ties can be physically understood quite naturally on the 
basis of a strong coupling expansion around the limit of 

Jl2 = |Jl3|=0. 

From the spectrum, we can also obtain the level- 
crossing fields at which the GS of the molecule changes 



from S =^toS 



and from S 



toS 



I We find: 



gtiBHci/kB ^ O.4OJ23, gtJiBHc2/kB ^ I.43J23. (9) 

With J23 = 217 K and the average g ~ 2.38 (see pre- 
viously), these numbers yield Hd ~ 54 T (compare to 
Fig. 14) and i?c2 -^ 194 T. 

GS properties. To probe the nature of the GS, we 
calculate several experimentally relevant GS expectation 
values. We first consider how the total S=\ moment of 
the GS (at low T, this moment can be attained by a small 
applied field) is distributed among the five Cu sites. At 
T<A, we find 



{St) = (Sf)o X tanh 



f 9lJ-BH \ 
\ 2kBT J ' 



(10) 



TABLE II. The exact spectrum of the spin Hamiltonian T-L 
of Eq. (8) Sit y — J12/J23 ~ f • For each state we also provide 
its total spin S, its parity p under the real-space inversion 
through the central site S3, as well as the expectation values 
of S22', S22', and S22/. The horizontal lines differentiate the 
three unperturbed manifolds of the strong coupling limit y — O 
(see text and Fig. 16). 
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FIG. 15. (Color online) Temperature dependence of the local 
moments of the three inequivalent Cu sites in CU5-NIPA as 
obtained by exact diagonalization using our ab initio exchange 
coupling parameters and a field of fioH = 1.6608 T. The 
inset shows the same results in a different scale, so that the 
deviation from the low-T behavior of Eq. (10), shown by the 
dashed lines, is better visible. Also, T* ~ 38 K denotes the 
characteristic temperature where the Cu2 moments change 
sign. 



where (<S'f)o are the GS expectation values (in an in- 
finitesimal field): 



(5|)o = 0.1416, (5|)o - (^|')o = -0-1415, 
(5f)o = (5f,)o = 0.3207. 



(11) 



The signs of the above numbers are in agreement with the 
classical ferrimagnetic GS discussed above (Fig. 1, right), 
whereby the spins S3, Si, and Si' are aligned parallel 
to each other, but antiparallel to S2 and S2' . However 
the strong deviation (especially for S3, S2 and S2') from 
their maximum (classical) value of i indicates that this strength of various spin-spin correlations e^ 



state has strong QM corrections, as can be directly seen 
from the explicit form of the wavefunction (not shown) . 
This feature will be explained on the basis of the strong 
coupling description developed below. 

The deviation of the local moments from Eq. (10) at 
high temperatures is shown in Fig. 15 at a field of 1.6608 
T. We note in particular the non-monotonic temperature 
dependence of the magnitude of the Cu2 moments, which 
was discussed in Sec. Ill C 1 above in relation to the NMR 
lineshape. 

To further probe the nature of the GS, we consider the 
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We find egi = egi/ ~ +0.2122, 632 = 632' ~ -0.4810, 
ei2 = 61-2' ^ -0.4299, ei2' = 61-2 ^ -0.2858, and 
fill' = 622' — +0.2345. These values corroborate the 
above QM ferrhxiagnetic picture for the GS, and in addi- 
tion highhght that the bonds with the largest exchange 
coupling J23 (i.e., the bonds 3-2 and 3-2') exhibit also 
the strongest correlations. Moreover, using the above 
correlations and the values of the exchange couplings, 
one finds that the J23-bonds contribute about 72% of 
the total GS energy Eq ~ -1.3289J23. 

We should remark here that the above result ew = 
622' does not arise from the symmetry of the Hamilto- 
nian, but it is an exact property of the particular GS 
only (see explanation below). 

Effective g-tensors. To make contact with our ESR 
findings in Sec. Ill B we calculate the effective g-tensors 
of the GS and the first excited state in terms of the in- 
dividual g-tensors of the three inequivalent Cu sites. To 
this end, one needs the equivalent spin operators within 
each multiplet in order to rewrite the Zeeman Hamilto- 
nian as 



^B-g, -S, K^B-geff-S 



(12) 



For the GS, the equivalent operators can be extracted 
from Eq. 11 above, which yields 

gcff, GS = +0.0708 g3 - 0.1415 g2 + 0.3207 gi . (13) 

So the deviation from the isotropic g — 2 value is coming 
mainly from that of gi . 

For the first excited state, our exact wavefunctions give 
the equivalent operators 
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leading to 



Global spectral structure. At y = 0, the spins Si and 
Si' are isolated and are, thus, free to point up or down. 
On the other hand, the spins S3, S2, and S2' form an 
AFM spin trimer with the Hamiltonian 



T~L{)/ J23 — S3 • (S2 



S2') = ^(SL'-SL-Sl)- (16) 



Hence the spectrum of 'Ho can be derived analytically in 
terms of the good quantum numbers S22' and S'322' : 



5*22', ■5'322', parity, deg, g^") 

1 \ even 8 — J23 

i odd 8 

1 I even 16 +J23/2 



(17) 



where we have also indicated the parity and the degen- 
eracy (deg) of each unperturbed level. The latter takes 
into account the four possible states of the space of Si 
and Si', and the Zeeman degeneracy. So we find three 
well separated unperturbed manifolds. The remaining 
couplings V = H — "Ho split these manifolds but, as we 
show below, most of these splittings are fairly weak, thus 
explaining the overall spectral structure at y = | (Ta- 
ble II). 

Since each unperturbed manifold has a well-defined 
S'322' , we can find the first-order splitting with the use 
of the "equivalent operators" : 
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(18) 



where the values of Ai^2,2' can be derived from the cou- 
pling (Clebsh-Gordan) coefficients of the respective an- 
gular momenta of each manifold (see below). Using 
A2 = A2', which holds in all manifolds, and replacing (18) 
in V. we obtain 



V ^ J.ff S, 



off ^3322' 



Sll' — TyJcR (S — S322 



'-S?i'), (19) 



off, 1st 



g3 + 0.3294 g2 + 0.004 gi 



(15) 



So here the deviation from the isotropic g = 2 value is 
mainly governed by that of g2 and g3. The contribution 
from gi is extremely small because in the first excited 
doublet the pair of Si and Si' forms almost a singlet 
(see Table II and Sec. IV C below). 

These expressions demonstrate why the effective g- 
tensors of the two ESR lines appearing below 70 K are 
different. 



C. Strong coupling description 



where Jeff = (A2 — ^3)y- Thus, to lowest order, the per- 
turbation gives rise to an effective exchange coupling be- 
tween Sii' and S322'. 

For the lowest eight-dimensional manifold, we couple 
the spin >S'3 = ^ to the spin S22' = 1 to get an S'322' = ^ 
object, which gives (see Appendix A): A3 = — ^, A2 = 
A2' =+|) find thus Jcff — y>0. The resulting correction 
to the lowest unperturbed manifold reads: 



S22', S322', Sii', S, parity p, deg, A^^^^ 



1 


1 

2 


1 


1 
2 


even 


2 


-y 


1 


1 
2 





1 
2 


odd 


2 





1 


1 
2 


1 


3 

2 


even 


4 


+2//2 



(20) 



We are now going to show that almost all properties 
of the model at j/ = J12 = —J13 = | can be adiabatically 
traced back - even on a quantitative level - to the strong 
coupling limit y = 0. 



Similarly, the second unperturbed manifold has S22' = 
and thus A2 = A2' =0 and A3 = l, giving Jcff = —y, which 
is now FM. So the first-order splitting of this manifold has 
the opposite structure from that of the lowest manifold. 
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FIG. 16. (Color online) Accuracy of the strong-coupling description. Evolution of the energy spectrum (a), the local moments 
(b), and the spin-spin correlations eij — {Si ■ Sj) (c) from the strong coupling limit y = up to y= |. The energies are given in 
units of J23. 



namely: 



5*22', 5'322', Sir, S, parity p, deg, A£:(^) 



odd 
even 
odd 





+y 



(21) 



The third manifold has ^322' — | 

A2' — 3, which in turn gives Jeff 

unperturbed manifold remains intact in lowest order, i.e., 

all 16 states have AE'^^) =0: 



and thus A3 = A2 = 
: 0. Thus, the third 



<S'22', <S'322', Sir, S, parity p, deg, AE'-^^ 



even 


2 





even 


2 





even 


4 





odd 


4 






(22) 



In particular, the energy of the state with the maximum 
spin (5 = |) remains equal to J23/2 in all orders of 
perturbation theory, as can be easily checked, e.g., for its 
maximum polarized M = | portion, which is a trivial 
eigenstate of H. 

Up to this lowest order, the level-crossing fields Hd 
and 7Jc2 a-re given by 

gUBHci/ks ^ 3y/2, gfJ.BHc2/kB ^ (3J23-y)/2. (23) 

With J23 = 217 K, y = 2J23/7, and the average g ~ 2.2, 
these numbers give: Hd — 54 T, and Hc2 — 170 T, which 
are close to the exact numbers at y = ^ given above. 

Figure 16 shows the evolution of the spectrum from 
y ~ up to y — ^, and demonstrates that the above 
first-order spectral structure reproduces quantitatively 
the spectrum at y = | . 

We should add here that 5*22', 'S'322', and ^n' do not 
remain good quantum numbers in higher orders of the 



perturbative expansion, which is expected, since the cor- 
responding angular momenta are not conserved quanti- 
ties under the full Hamiltonian. Still, their expectation 
values reported in Table II are very close to the above 
strong-coupling values. 

As a final remark, we note that one may use the above 
basis provided by the lowest-order theory and the good 
quantum numbers S, M, and p to find the invariant 
subspaces of H. Most of the subspaces turn out to be 
two-dimensional, except for the space with S = M = ^ 
and p =even (in which the GS belongs), which is three- 
dimensional. Here, it may be readily checked that each 
of the three states has 5*11' = 5*22', thus explaining the 
equality en' = 622' discussed above. The excited states 
do not share this property, so at higher temperatures the 
two correlations deviate from each other. 

GS properties. Let us now return to the nature of the 
GS. According to the above lowest order result, the GS is 
a total S — ^ doublet with even parity and total energy 
E'-^^ — — J23 — J12 • To understand the nature of this state, 
we note that it arises from the coupling of a spin s = 5 
(here S322') and a spin s = 1 (here Sn') to a total S—^ 
state. So we may again use the "equivalent operators" 



'322' 



-S, Si, Si' — > + — S, 



(24) 



which in conjunction with the ones relating Si, S3, and 
S5 to S322' (see Eq. (18) above) gives the central rela- 
tions: 



S2, S2 



_2o 

9=7 



Si, Si 



-is. 



(25) 



This equation provides a very clear explanation for 
the peculiar local moment distribution discussed above, 
since, when the total S — h moment is saturated (by 
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applying a small field), this equation gives: 



(SI 



(Si,) = -i, (Sf) = (S^,) 



(26) 



The same numbers can be extracted directly from the 
explicit form of the GS wavefunction, which for the M ~ 
+ ^ portion reads: 

lt)=-^lt)3|tl)ll'|t-l)22' + ^|t)3|<0)ll'IM22' 

+ ^i;)3(|ii)ii'IM22' - |io)ll'|tl)22'), (27) 

where \ti)ij = Ititj), \t-i)ij = \iilj), and |io)jj = 
(lUj) + \Uj)) /%/2- The leading term of Eq. (27)) is 
the classical ferrimagnetic configuration. However, the 
overall form of the wavefunction demonstrates explicitly 
the presence of strong QM corrections to the GS. 

The strong coupling limit reproduces the ferrimag- 
netic arrangement of the local moments and even their 
strengths, except for S3 which has almost tripled its 
moment at y = i. This is further demonstrated in 
Fig. 16(b), which shows the evolution of the local mo- 
ments from y = up to y= |. All moments, except (S*!), 
exhibit a small slope, i.e., higher-order corrections are 
weak. 

In a similar fashion, the strong coupling prediction 
for the spin-spin correlations are: 631 — e^y — g, 632 — 
632' = -5, ei2 = ei'2' = -5, ei2' = ei'2 = -|, and 
eii' = 622' = J I which are in close agreement with the 
exact values for y = f given earlier. This is again demon- 
strated in Fig. 16(c) that shows the evolution of the spin- 
spin correlation strengths from y = up to y = | . All 
correlations exhibit weak corrections, except ei2 which 
shows a somewhat larger slope in y. 



V. DISCUSSION AND SUMMARY 

The CU5-NIPA molecule comprises two non-frustrated 
spin triangles with two AFM couplings and one FM 
coupling each (Fig. 1, left). This coupling regime is 
rather unusual and deserves a further analysis with re- 
spect to the underlying crystal structure and interacting 
orbitals. Three leading couplings - Ji2,Ji3, and J23 - 
represent the Cu-O-Cu superexchange running via the 
//3-O atom in the middle of the Cua triangle. According 
to Goodenough-Kanamori rules, high values of the bridg- 
ing angle should lead to an AFM interaction, whereas low 
bridging angles close to 90° favor FM couplings. Our re- 
sults for CU5-NIPA (Table I) follow this general trend. 
However, a closer examination pinpoints two peculiari- 
ties of this compound. 

First, the coupling remains FM for the bridging angle 
of 107.9°, although the standard and commonly accepted 
threshold value of the FM-AFM crossover is slightly be- 
low 100°.^" Second, Jf^ is much larger than J™, even 



though the respective bridging angle is also larger and 
should lead to a smaller FM contribution. Both pecu- 
liarities should be traced back to the twisted configu- 
ration of the interacting CUO4 plaquettes (Fig. 1, left). 
The systematic work on the angular dependence of the 
exchange coupling is usually restricted to systems with 
two plaquettes lying in the same plane or only weakly 
twisted (the dihedral angle between the planes is close to 
180 deg). The Cu5 molecule represents an opposite limit 
of strongly twisted CUO4 plaquettes, with the dihedral 
angles of 91.8 deg (J12), 125.2 deg (J13), and 99.9 deg 
( J23) between the CUO4 planes. ^^ 

The large FM contribution to J12 (Ji™ = -134 K) 
can be ascribed to the nearly orthogonal configuration 
of the interacting CUO4 plaquettes. However, the mech- 
anism of this FM interaction is yet to be determined. 
The large FM coupling for the bridging angle of 90 deg 
is generally understood as the Hund's coupling on the 
oxygen site^^ or the direct FM exchange between Cu and 
O.^'' Both mechanisms depend solely on the Cu-O-Cu 
angle and should be rather insensitive to the mutual ori- 
entation of the CUO4 plaquettes. Therefore, other ef- 
fects, such as the direct Cu-Cu exchange and the Cu- 
O-Cu interaction involving fully filled Cu 3d orbitals, 
may be operative. Experimental information on the mag- 
netic exchange between the strongly twisted CUO4 pla- 
quettes remains scarce and somewhat unsystematic. The 
Cu-O-Cu angle of 104.5 deg combined with the sizable 
twisting give rise to a FM nearest-neighbor interaction 
in the kagome material kapellasite Cu3Zn(OH)gCl2,^^ 
whereas a similar geometry with the Cu-O-Cu angle of 
107.6 deg results in an overall AFM coupling in dioptase, 
CueSieOis • 6H20.^5 

We also note that CU5-NIPA brings a fresh perspec- 
tive on magnetostructural correlations in Cus triangular 
molecules, where the Cu-X-Cu bridging angle was previ- 
ously considered the key geometrical parameter. As long 
as the ligand atom X lies in the CU3 plane, the bridg- 
ing angles remain close to 120 deg and should generally 
lead to the AFM exchange. The FM exchange would 
only be possible when the ligand atoms are shifted out of 
the plane, as in the CI- and Br-containing Cus molecules 
where the Cu-X-Cu angles are below 90 deg.^® The twist- 
ing of the CuX4 plaquettes provides another opportunity 
for creating the FM exchange and, moreover, for intro- 
ducing it selectively. Surprisingly, organic ligands have 
only a weak effect on the magnetic interactions. Al- 
though carboxyl groups (COO^) provide an additional 
superexchange pathway between Cul and Cu3 (Fig. 1, 
left), their molecular orbitals do not influence the Cu 
dx2_y2 Wannier functions (Fig. 12). For example, the 
FM contribution to J12 exceeds that to J13. This demon- 
strates that the FM exchange is basically unrelated to the 
organic ligands. It should be understood as a joint effect 
of the the low Cu-O-Cu angle and twisting. 

The combination of FM and AFM couplings also has 
an important effect on the magnetic GS of CU5-NIPA. 
A regular spin triangle entails the four-fold degenerate 
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GS that further sphts into two close-lying doublets by 
virtue of anisotropy,^^ residual interactions between the 
triangles/^ or a marginal distortion of the triangle.^* In 
CU5-NIPA, there is only one doublet state, which is sep- 
arated by about 68 K from the first excited state (see 
Table II). According to Kramer's theorem, in zero field 
the degeneracy of this state can not be lifted, hence the 
tunnel splitting is exactly zero, and no magnetization 
steps (Landau-Zener-Stiickelberg transitions) should oc- 
cur, in contrast to, e.g., the Ve molecule. ^° Therefore, 
broad butterfly hysteresis efi^ects but without tunneling 
are expected in the magnetization process of CU5-NIPA. 

Of particular interest is our experimental finding of the 
enhanced ^H nuclear spin-lattice relaxation rate 1/Ti at 
a characteristic temperature slightly below the spin gap 
(T ~ 40 K). While such an enhancement has been re- 
peatedly found in numerous AFM homometallic'^^ and 
heterometallic'^'^ rings of spins S > |, it is very rare for 
S = -^ and has been observed only in molecules with a 
high-spin GS.^"^'®^ The origin of the peak is likely the 
same in both cases, namely, the slowing down of the 
phonon-driven magnetization fiuctuations.'^*''*'^"''^ How- 
ever, the sparse excitation spectrum and the presence 
of nonequivalent Cu sites with different local magnetiza- 
tion render CU5-NIPA dissimilar to typical homometallic 
rings. We argue that in CU5-NIPA inter-multiplet Or- 
bach processes make the dominant contribution to the 
spin-lattice relaxation process, as in the heterometallic 
ring CrrNi at very low temperatures.^^ 

The use of molecular magnets in quantum computing 
is severely restricted by their short coherence time. Co- 
herence times on the order of 100 ^s could be achieved 
in, e.g., V15 by arranging magnetic molecules in a self- 
assembled layer formed by an organic surfactant. ^^ This 
method is fundamentally similar to the formation of the 
metal-organic framework in CU5-NIPA. Therefore, it may 
be interesting to study the decoherence process in this 
molecular magnet and, more generally, explore the role 
of organic bridges in the decoherence process. 

Finally, we would like to note that only a few examples 
of magnetic pentamers based on spin-^ ions have been re- 
ported so far. The [Cu(/i-L)3]2Cu3(Ai3''-OH)(PF6)3-5H20 
{L~ = 3, 5-bis(2-pyridil)pyrazolate) shows a somewhat 
similar phenomenology with the spin S = ^ GS and the 
respective magnetization plateau ranging up to 30 T.^^ 
However, the geometrical structure of this compound is 
a trigonal bipyramid that is notably different from the 
hourglass shape of the magnetic cluster in CU5-NIPA. 
The magnetic pentamer based on Cu5(OH)4(H20)2(A- 
a-SiW9033)2] is more similar to our case and also reveals 
the spin S — ^ GS.^^'^* Nevertheless, the equivalence of 
J23 and J13, as well as sizable long-range couplings, dis- 
tinguish its spectrum and magnetic properties from that 

of CU5-NIPA. 

In summary, we have studied thermodynamic proper- 
ties, spin dynamics, microscopic magnetic model, and en- 
ergy spectrum of the spin-^ Cus pentamer in GU5-NIPA. 
This magnetic molecule has an hourglass shape with two 



AFM couplings and one FM coupling on each of the two 
triangles that form the non-frustrated pentamer. In zero 
field, the ground state of GU5-NIPA is a doublet with the 
total spin oi S = | . This ground state is separated from 
the first excited state by an energy gap of A ~ 68 K. Our 
results evidence a highly inhomogeneous distribution of 
magnetization over rhe nonequivalent Cu sites accord- 
ing to the quantum nature of the magnetic ground state. 
The maximum in the spin-lattice relaxation rate 1/Ti is 
very rare among molecular magnets with spin- ^ ions and 
a low-spin ground state. 
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Appendix A: Coupling of a spin 1 to a spin 1/2 
object and "equivalent operators" 

Here we work out the coupling of two spins, one with 
spin Sa — 1/2 and the other one with She — 1- For the 
latter, we shall imagine that there are two spins Si, = 
Sc = 1/2 forming a triplet, namely |ii) = |tt), l^-i) = 
III), and \to) = -Tjdtlj + llt))- The states of the system 
can be labeled by \S, M), where S is the total spin (here 
S = i or |), and M is the projection along some axis z. 
The states with S = ^ are given by: 



i> = |t)a®|tl)6c, 



'2' 2^ V3 



-^ {\ i) a <»\tl} be + V2\t) a <»\ta)b. 

\t)a <» \t^l)bc + V2\i)a <» \to}b. 



-1/fcc- 



On the other hand, the two states of the S — ^ doublet 
read: 



|i i\ 
I2' 2/ 



1 

7! 



1 



V3 



\t) a <» \to) be -V2\i} a <»\tl} be), 
-\i)a<»\to)be + V2\t)a®\t-l)be 



Using the above explicit relations, it is straightforward to 
find out the equivalent operators within the 2x2 manifold 



of the 5' = i doublet. We find: 
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Sa 



rS, Sf, 



, ^J ^c 



(Al) 



Appendix B: Structural data 

Table III shows the equilibrium positions of the 17 hy- 
drogen nuclei in the structure. 



TABLE III. Relaxed hydrogen positions in Cus-NIPA. Last 
column lists atoms bonded to the given hydrogen atom. 
The notation of atoms follows Ref. 22, a = 10.8303 A, 
h = 11.4692 A, and c = 11.5697 A. 





x/a 


y/h 


z/c 




H2 


0.0341 


0.9481 


0.3828 


(C2) 


H4 


0.4295 


0.6383 


0.4380 


(C4) 


H6 


0.4442 


0.0066 


0.2381 


(C6) 


HIO 


0.4563 


0.4564 


0.2062 


(CIO) 


H12 


0.1278 


0.8023 


0.0923 


(C12) 


HIS 


0.0814 


0.3482 


0.1695 


(013) 


H14 


0.5705 


0.7834 


0.9955 


(C14) 


H14A 


0.1105 


0.1813 


0.9704 


(014) 


H14B 


0.0012 


0.2993 


0.9133 


(014) 


H15A 


0.6634 


0.2163 


0.2633 


(015) 


H15B 


0.8008 


0.0961 


0.2947 


(015) 


H16A 


0.9178 


0.6522 


0.4705 


(016) 


H16B 


0.8144 


0.5794 


0.5639 


(016) 


H17A 


0.3996 


0.2783 


0.2640 


(017) 


H17B 


0.4504 


0.1667 


0.3606 


(017) 


H18A 


0.0954 


0.1403 


0.6165 


(018) 


H18B 


0.1939 


0.1190 


0.6921 


(018) 
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